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ABSTRACT 

Wc present a comprehensive spectroscopic imaging survey of the distribution 
and kinematics of atomic hydrogen (HI) in 16 nearby spiral galaxies hosting low 
luminosity AGN, observed with high spectral and spatial resolution (resolution: 
~20", ~5 km s~^) using the NRAO Very Large Array (VLA). The sample con- 
tains a range of nuclear types, ranging from Seyfert to star-forming nuclei and 
was originally selected for the NUclei of GAlaxies project (NUGA) - a spectrally 
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and spatially resolved interferometric survey of gas dynamics in nearby galaxies 
designed to identify the fueling mechanisms of AGN and the relation to host 
galaxy evolution. Here we investigate the relationship between the HI proper- 
ties of these galaxies, their environment, their stellar distribution and their AGN 
type. The large-scale HI morphology of each galaxy is classified as ringed, spiral, 
or centrally concentrated; comparison of the resulting morphological classifica- 
tion with AGN type reveals that ring structures are significantly more common in 
LINER than in Seyfert host galaxies, suggesting a time evolution of the AGN ac- 
tivity together with the redistribution of the neutral gas. Dynamically disturbed 
HI disks are also more prevalent in LINER host galaxies than in Seyfert host 
galaxies. While several galaxies are surrounded by companions (some with asso- 
ciated HI emission), there is no correlation between the presence of companions 
and the AGN type (Seyfert/LINER). 

Subject headings: galaxies: active - galaxies:kinematics and dynamics - galax- 
ies:Seyfert-galaxies:individual:(NGC1961, NGC2782, NGC3147, NGC3368, NGC3627, 
NGC3718, NGC4321, NGC4569, NGC4579, NGC4736, NGC4826, NGC5248, 
NGC5953, NGC6574, NGC6951, NGC7217)- radio hues: galaxies 



1. Introduction 



Active Galactic Nuclei (AGN) represent some of the most extreme conditions and include 
the most powerful individual objects to be found throughout the Universe. Nowadays, the 
phenomena of nuclear activity is generally understood to be the result of accretion of material 
onto a SuperMassive Black Hole (SMBH); very likely through the infall of gas from its host 
galaxy (e.g. Rees||1984 ). Although SMBHs exist in most galaxies, only a small fraction of all 
galaxies exhibit nuclear activity ( Huchra fc Burg|1992 Ho et al. 1997 Miller et al.|2003 ). The 
majority of AGNs show low-ionization narrow emission-line regions (LINERs), which have 
luminosities lower than Seyfert galaxies and quasars. Originally, LINERs and Seyferts were 
considered to form a continuous distribution in the usual line-ratio classification diagrams. 
Recent investigations of the host properties of emission-line galaxies from the Sloan Digital 



Sky Survey (Kewley et al. 2006) confirmed that their nuclear activities are due to accretion 



by a SMBH. But interestingly, LINERs and Seyferts are also clearly separable in emission line 



ratio diagrams (Groves et al. 2006), and that these two classes have distinct host properties 



( Groves et al.|200"6 ), implying that their distinction is in fact not simply an arbitrary division. 
An open question is therefore what mechanisms give rise to the different AGN types observed. 



A coherent picture of AGN feeding, e.g. how to remove the angular momentum from 
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host galaxy gas so that it can reach the central parsec, is still missing. Therefore it is ex- 
pected that a hierarchy of mechanisms might combine to transport gas from large kpc scales 
down to the inner pc scales dShlosman et al.||1990t|Combesl2003||Jogee ||2004D . Additionally, 
trigger mechanisms driven by the galactic environment, like minor mergers and tidal forces. 



might play a role as well (Barnes & Hernquist 1991). A larger gas fraction for Seyferts than 



for normal galaxies was suggested by Hunt et al. ( 1999 ) and a prevalence of stellar rings (RC3 
classification) in active galaxies (Seyfert and LINERs) compared to non-active galaxies was 



suggested (Hunt & Malkan 1999). NICMOS imaging of the centers of 250 nearby galaxies 
( Hunt fc Malkan||2004 ) revealed systematic morphological differences between Hll/starburst 
(most disturbed), Seyfert (intermediate disturbed), LINER and normal galaxies (most reg- 
ular). The study of the circumnuclear regions of 24 Seyfert2 galaxies using HST images 



(Martini & Pogge 1999) revealed a dominance of nuclear spirals, suggesting that spiral dust 



lanes may be responsible for feeding gas to the central engines. 

Theoretical simulations have also made progress in addressing the nature of nuclear 
fueling with respect to different types of gravitational instabilities and their feeding efficiency: 
nested bars (e.g. Shlosman et al.|[T989 Friedli fc Martinet||1993 Maciejewski fc Sparke||2000 



Englmaier & Shlosman'2004j), gaseous spiral density waves (e.g. Englmaier fc Shlosman|2000 



Maciejewski||2002 2004aj|bj), m = 1 perturbations (e.g. Shu et alT]|1990 Junqueira & Combes 



1996 Garci'a-Burillo et al. 2000) and nuclear warps (Schinnerer et al. 2000) have all been 



suggested as possible transport mechanisms. 

In order to distinguish models for nuclear fueling, observations of neutral gas with high 
angular and velocity resolution are required. Therefore, the IRAM key project NUclei of 
GAlaxies ( NUGA; see Garcia-Burillo et al. 2003 ) was established - a spectroscopic imaging 



survey of gas in the centers of nearby low luminosity AGN. As most of the gas in galaxy nuclei 
is in the molecular phase, the survey used millimeter CO lines to conduct a detailed mapping 
of molecular gas dynamics at high-resolution (0.5") in the central kiloparsec of AGN hosts. 
The CO-NUGA survey (using the IRAM Plateau de Bure Interferometer) reveals a wealth 
of nuclear gas distribution and kinematics which are studied in detail: (a) m=l gravitational 



instabilities (one-arm spirals and lopsided disks; NGC 4826: Garcia-Burillo et al. (2003)) 



(b) m=2 instabilities (two-arm spiral wave) expected to form in stellar bar potentials show 



only small amounts of molecular gas coincident with the AGN (e.g. NGC 4569: Boone et al. 
( |2057| , NGC 4579: |Garcfa-Burillo et al.| ( |2005j ) NGC 6951: Schinnerer et al., in prep.), (c) 
stochastic spirals that are related to non self-gravitating instabilities, and rings (NGC 7217: 



Combes et al. 


2004 


), and (d) lar 


(NGC 3718: 


Krips et al. 2005D 



components resembling jet emission from AGN were detected (Krips et al. 2007a) and a 



molecular gas disk/torus of dense gas was observed in the HCN line emission (NGC 6951: 
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Krips et al. 2007b). Additionally, Garci'a-Burillo et al. (2005) derived in a pilot study the 



gas inflow rates via measurements of the gravitational torque onto the molecular gas in the 
central kiloparsec for 4 galaxies. 

To complement these CO data and provide a more complete census of the ISM and 
gas flows from the outskirts of the galaxy disks to the very center, the HI-NUGA project 
was initiated. HI-NUGA provides observations of the HI gas distributions and kinematics 
for 16 galaxies (including archival data for 2 galaxies) of the NUGA sample using NRAO's 
Very Large Array (VLA). Determining the HI distribution and kinematics is necessary to 
establish whether characteristics of the nuclear dynamics, and hence the nuclear activity, 
depend on overall host galaxy properties. Possible dependencies might also exist between 
nuclear modes and large scale drivers (i.e. tidal forces exerted by companions), or the overall 
content /distribution of the atomic gas. Since the HI gas extends in most disk galaxies much 
further than the optical disk, it is more loosely bound in the outer region because of the 
decrease of the gravitational potential. Hence, it provides an ideal tool to identify interactions 
and tidal features (e.g., Simkin et al. 1987 Mundell et al. 1995). Furthermore, due to its 



dissipative nature the gas is more sensitive to dynamical disturbances, both internal ones 
such as non-axisymmetric potentials or external ones (e.g. tidal interactions). Thus, one 
could expect, for example, to find a prevalence of Seyfert activity with asymmetries in the 
gas distribution and kinematics. 

In this paper we present for 16 nearby active galaxies data obtained in the 21 cm 
emission line of neutral hydrogen using the VLA. The relation between HI properties, such 
as the large scale environment (Hi/optical companions, HI disk asymmetries) and the AGN 
type (Seyfert, LINER, HII) are analyzed in detail. The outline of the paper is as follows: 
We describe the HI-NUGA sample, the observation and data reduction in ^ The results 
and derived HI properties are presented in ^ along with a discussion of the presence of 
companions and tidal disturbances of the HI disks and any correlation with AGN type. 
A comparison of AGN types and optical (stellar) light distribution is also examined. The 
results are discussed in the context of correlations with the HI environment (e.g. tidal forces) 
and the AGN type (Q. A summary is presented in ^ 



2. Sample Description, Observations and Data Reduction 
2.1. Sample Selection and Survey Goals 

The NUGA project is a multiwavelength study of the nuclei of nearby AGN. The core 



sample ( Garcia-Burillo et al. 2003) consists of 12 nearby AGN with D < 54Mpc selected 
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for the availability of high-quality optical and near infrared (NIR) images obtained with 
both ground-based telescopes and the Hubble Space Telescope (HST). The distribution and 
kinematics of molecular gas in all galaxies in the sample have been observed with the IRAM 
PdBI at 1mm and 3mm with maximum angular and spectral resolutions of 0.5" and 3.6 km 
s~^, respectively. The NUGA supersample adds another 16 galaxies with CO data of similar 
quality obtained by the NUGA team with the IRAM array as part of other projects. 

For our detailed HI study using the VLA, 16 galaxies were selected (15 of which have 
high quality CO data of < 1" resolution): the entire CO-NUGA parent sample, 3 additional 
galaxies from the CO-NUGA supersample and NGC 5248. The HI-NUGA sample therefore 
covers a variety of nuclear activity: 7 Seyferts, 7 LINERs and 2 HII galaxies. All galaxies 
in our sample are spiral galaxies, ranging in Hubble type from Sa to Sbc in disk type. 
Typical spatial resolution of the HI data is 21"=(1.9 kpc) for a mean distance of ~20 Mpc; 
the distance of our sample ranges from 4 Mpc to 54 Mpc. An overview of the general 
properties of our HI sample is given in Tab. [Tj The overall aims of the HI-NUGA project 
are to 1) determine the global HI properties of the sample galaxies and their large scale 
environment to search for relations between the nuclear activity and the HI environment, 
for all disk properties of the full sample, 2) to trace the gas inflow on various spatial scales 
for the best candidates of the HI-NUGA sample with high angular resolution HI data (VLA 
B-I-CD array) and to determine the fueling rates by using gravity torques, 3) to obtain 
detailed dynamical models (using SPH codes with a fixed gravitational potential and/or 
N-body codes.) describing the gas flow in the galactic disk over several radii for the most 
representative galaxies. In this paper we concentrate on goal 1). 



2.2. AGN classification 



Active galaxies can be separated into the categories of emission-line AGNs and HII 
regionlike galaxies (starburst) based on the differing nature of the photoionization source. 
Emission-line AGNs can be divided into Seyferts and LINERs. The majority of AGN are 
LINERs with nuclear optical spectra that are dominated by emission lines of low-ionization 



species such as [01] A6300, [Oil] A3726,9 and [SII] A6717 ( |Ho et al.|1997| ). The main difference 
between Seyferts and LINERs is thought to be accretion rate, which explains why LINERs 
have lower AGN luminosities than Seyfert galaxies. Although Seyferts and LINERs are 
basically powered by the same engine, i.e. accretion onto a SMBH, the Seyfert/LINER 



distribution is bimodal (Groves et al. 2006), implying different accretion mechanisms (for a 



Ulrich et al. 


1997 


Combes 


2001 



LINERs seem in general to be hosted by older, more massive, less dusty galaxies with 
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higher stellar velocity dispersions and lower nuclear [OIII] luminosities than Seyfert galaxies 
(Kewley et al. 2006). Despite numerous studies in this field of research, the origin of their 



differences (e.g. different stages in evolution) is still under debate. 

Although all the galaxies in our sample have several nuclear classifications available in 
the literature, we reclassified their nuclear types as described in the following, in order to have 
a coherent scheme for our sample. The most common method to select star-forming galaxies 



is based on the Baldwin, Phillips & Terlevich (1981) empirical diagnostic diagrams using the 



optical line ratios of [OI]/Ha, [SII]/Ha, [NII]/Ha, and [0III]/H/5. For the purpose of this 
paper we separated our sample in Seyferts, LINERs and non-AGN (HII+Transitions) objects 
following the most recent method of Kewley et al. (2006): First we distinguished AGN 



galaxies from non-AGN galaxies by using the optical line ratio [0111] /H/5 versus [Nil] /Ha 
(left panel Fig. [T|. The line ratios are taken from the cataloged data of Ho et al. (1997) 



and Kim et al. (1995). For NGC 6574 no accurate line ratios were available in the literature 



for this method and we used its classification listed in NED as Seyfert. The line separating 



AGN and Starburst/Transition galaxies is given by (Kewley et al. 2006): 

[OIII]. 0.61 



log{- 



Ha 



-0.47 



1.19 



(1) 



Secondly, we divided the AGN galaxies in Seyfert and LINERs by using the [0III]/H/9 versus 
[01] /Ha diagram (right panel Fig. [T]) with a separating line defined by: 



log{ 



(2) 



The classification of our sample results in 2 non-AGN galaxies (HII and transition objects), 
7 Seyferts and 7 LINERs. Since only 2 non-AGN galaxies are present in our sample, we will 
mainly focus on a comparison between the properties of Seyferts and LINERs. A detailed 
comparison of HI properties between Seyfert and non-AGN galaxies is carried out by Mundell 
et al. (in prep.). The average morphology type taken from Hyperleda is roughly the same 
for Seyferts (t=2.81) and LINERs (t=2.69). Hence, effects due to bias of the morphology 
type are unlikely. 



2.3. VLA Observations and Data Reduction 

The 16 galaxies presented here were observed in 2003 and 2004 in their HI line emis- 
sion using the VLA in its C- and D-configuration. Three galaxies (NGC 3627, NGC 4736, 
NGC 4826) were also observed as part of the THINGS project (Walter et al., in prep., BC- 
array data) while adequate data for NGC 2782 and NGC 3718 were already available in the 
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VLA archive. For NGC 3718, only the C-array data has been analyzed as the D-array data 
was lacking good phase calibrator measurement required for our purposes. 

Flux calibrator measurements were performed at the beginning and at the end of each 
observation cycle. The phase calibrator was observed before and after each source cycle with 
a maximum distance between source and phase calibrator of 12°. The correlator spectral 
setup used was spectral line mode 4 with Hanning smoothing and 64 channels per 1.5625 MHz 
channel width per frequency band providing a frequency resolution of 24.414 kHz/channel 
(~ 5.2 km s"^). NGC 2782 was observed with a resolution of only 48.83 kHz (~ 10.48 km 
s~^). Our observations have an average on-source integration time of 2.7h in the C-array 
and 2.8h in the D-array configuration. Thus, the corresponding RMS flux sensitivity of 0.6 
and 0.8 mJy /beam/channel gives a 3-(Tdetection limit of ~ 4 x 10^^ cm~^ and ~ 0.38 x 10^^ 
cm~^ column density for the C- and D-array, respectively. 

The data reduction was performed following the standard procedure of the Astronomical 



Image Processing System (AIPS; Greisen 2003). Calibration solutions were first derived for 



the continuum data-set (inner 3/4 of the spectral band width) and then transferred to the 
line data. The bandpass solutions were determined from the phase calibrator measurements 
to account for channel to channel variations in phase and amplitude. 

After the calibration was applied, the two-IF data set was split into single-source data 
sets with 64 channels each. For each set the first 8 and last 8 channels were discarded, 
resulting in two 48-channel data sets. These were then combined into a final cube of 96 
channels with an effective bandwith of 2343.7 kHz (~ 500 km s~^) and a channel width of 
24.4kHz (~ 5.2 km s'^). For the THINGS data and NGC 2782 only the first 4 and last 
4 channels were discarded which resulted in 112-channel sets with an effective bandwith 
of 2734.4 kHz and 5469kHz (NGC 2782), respectively. For NGC 1961 two independent 
frequency setups were used with an effective bandwith of 4687.5 kHz (~ 1018 km s~^). After 
that, the cahbrated C and D array data cubes were combined. 

The task UVLIN within AIPS was used to subtract the continuum emission. The 
continuum emission was measured only in channels without line signal (usually the first 10 
and last 10 channels of a data cube). The fit was done with a first order polynomial, namely 
a DC term (offset from zero) and a slope. For NGC 5953 only the first 8 channels were used 
(simple offset fit) because of additional line signal in the last channels. Channels affected by 
galactic HI emission (NGC 4569, Channel 4-7) and internal RFI (NGC 1961, Channel 17-18) 
were dismissed and have been excluded from the further analysis. However, these channels 
were well outside line emission in both cases. 



After the continuum substraction, the data were CLEANed using the task IMAGR 



- 8 - 



within AIPS. The cahbration of each observation (C and D configuration) was verified indi- 
vidually (by a check of the imaged cubes). The fiux-limit for CLEANing was set to a peak 
value in the residual image of 2xRMS with a maximum iteration number of 10^ steps. After 
CLEANing the residuals were rescaled by the ratio of the restoring beam area to the dirty 
beam area. The reason for this rescaling is that typically the residual fiux of the source in 
cleaned channel maps is overestimated due to the different shape of the dirty and CLEAN 
beams. However, the results of the rescaling method turned out to not be significantly differ- 
ent from the results obtained without rescaling (less than 5%). We produced two CLEANed 
data cubes for each galaxy: 1.) Naturally weighted imaging with a velocity resolution of 
~20.8 km s~^ and an average angular resolution of 33", in order to obtain a datacube with 
maximum sensitivity to the outer HI disk and for HI companions around it. 2.) Robust 
weighted imaging for a high angular (20") and spectral (~5.2 km s~^) resolution data cube 
to study the galactic HI disk itself. To find the best compromise between angular resolution 
and RMS, several robust weighting parameters were tested and a robust weighting parameter 
of 0.3 was selected. The cellsize of each grid was set to 4"/pixel for natural and 2"/pixel 
for robust weighting with a resulting mapsize of 1024 x 1024 pixel. That corresponds to a 
field of view of ~ 68' (natural) and ~ 34' (robust) such that the 30' primary beam is fully 
mapped in both cases. The RMS values and beam sizes for natural/robust weighting are 
listed in Table [2] with an average achieved RMS value of 0.34/0.46 Jy beam~^ respectively. 
The channel maps for the naturally weighted cubes are presented in the appendix (see Fig. |2] 
as an example). All further analysis was performed on these data cubes. A primary beam 
correction was only applied for the fiux measurement. 

The subsequent analysis has been done with the Groningen Image Processing SYstem 



(GIPSY; [van der Hulst et al. 1992). The channel maps were combined to produce zeroth 
(intensity map), first (velocity field) and second (dispersion field) moments of the line profiles 
using the task MOMENT. The RMS values have been measured in two regions where no 
line signal was apparent and averaged over all channels of the cubes. A fiux cut-off of three 
times the channel-averaged RMS value was used for the moment maps. In order to exclude 
one channel wide noise peaks from the moment maps, two adjacent channels with fiux values 
above the cut-off level were required (corresponding to a width of 10.4 km s~^ for the robust 
weighted cubes). However, for the natural weighted data cube we did not apply such a 
requirement since a single channel has already a width of 20.8 km s"^, which is also the 
range of the minimal velocity width of an expected HI emission line. 
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3. Analysis and Results 

In the following the distribution and kinematic properties of the HI gas disks of the 
Hl-NUGA sample will be derived, compared to optical properties of the host galaxies and 
analyzed as a function of the AGN type. 



3.1. HI Distribution and Environment 

3.1.1. HI Gas Morphology 

The velocity-integrated HI intensity maps and the optical light distribution (DSS2-blue 
maps) are presented in Fig. |3] The HI intensity maps show a large variety of morphologies 
in the disk gas distribution including spiral (e.g. NGC 6951), ringed (e.g. NGC 7217) and 
centrally concentrated geometries (e.g. NGC 5953). A simple classification by eye has been 
used to assign a ring or spiral morphology. In cases where the distinction was unclear, 
the intensity map was projected into polar coordinates, giving the azimuthal variation of 
the intensity around a given center. If the intensity was constant along a radius, a ringed 
geometry was assumed. An overview of the Hl-morphology of our sample is presented in 
Tab. [3j In summary, we find spiral geometry in 8 galaxies, ringed geometry in 4 galaxies, 
and concentrated geometry in 5 galaxies, whereby NGC 4826 shows concentrated/spiral 
geometry and NGC 3368 ringed/spiral geometry in the central/outer disk. For NGC 3718 
no classification is possible as its disk is nearly edge-on and shows in addition a large warped 
disk. 

In the following, we describe the HI morphology of the disk for the three galaxies which 
clearly differ from the rest of the HI-NUGA sample: NGC 5248 is the only galaxy of our 
sample which shows a large-scale HI "gas bar" (13 kpc in length). This HI bar coincides with 
a stellar bar (see left panel of Fig. [3]). In addition, NGC 5248 shows two large HI spiral arms 
located outside the HI bar itself. Thus its HI morphology is remarkably similar to NGC 4151 



which has a gas-rich bar and optically faint but gas-rich outer spiral arms (Mundell et al. 



1999). 



NGC 3368 exhibits two rings with radii of about 78 " (~3.1 kpc) and 220 " (~8.6 kpc) 
for the inner and outer ring, respectively. The outer ring could also be two closely wound 
spiral arms which are not resolved at our resolution of 19". The inner ring is ~15 times larger 



than the small-scale ( r ^ 5") nuclear "mini-bar" (seen in the optical and NIR, Moiseev et 



al. 2004), which extends along a position angle of 120°— 125°. 



NGC 4826 has two counter-rotating HI gas disks ( Braun et al.||1992 ) where the inner HI 
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disk has a radius of ~70 " (~1.4 kpc) and the outer HI disk a radius of ~700 " (~13.9 kpc). 
No obvious asymmetries are seen, neither in the inner disk nor in the outer disk. Here the 
radius of the optical disk (~300") hes between the radii of the two counter-rotating parts of 
the HI disk. 

Two galaxies in our sample show very distorted outer HI distributions in the form of HI 
tails and large faint extensions stretching from only one side of the main galaxy. NGC 2782 
has a large HI tail with a projected length of ~ 6' (~3 times the HI radius, see 3.2 and 



Smith 1991, 1994) extending from its center to the north-west. NGC 1961 shows a large HI 
extension to the north-west with a projected length of ~ 4' (~ 2 times the HI radius; see also 



Shostak et al. 1982) from the center. NGC 3718 exhibits a large HI warp (Schwarz 1985) 



Finally, two galaxies show large differences between the distribution of the stellar light 
and the HI gas: For NGC 6574 the HI gas distribution shows two HI peaks inside the stellar 
disk well separated by 30" even at our angular resolution of 19" (see Fig. |3]). NGC5953 
forms a close binary (likely merging) pair with NGC5943 (see left panel Fig 3) with the two 
galaxy centers separated by only 60" and a difference in systemic velocity of ~53 km s^^. 
The HI distribution shows a wing-like structure whose velocity field could be interpreted as a 
single disk, roughly uniformly distributed with a faint HI extension to the north. Also the HI 
velocity field shows no obvious tidal asymmetries which could be associated with the merging 



of the two optical galaxies. Our study of the mass content (described in ^3.4) revealed a 
larger value for NGC 5953 {Log{MHj/ D"^) = 7.1) than the typical HI mass content (Bettoni 



et al. 2003) for a galaxy of that Hubble type (Log(M////D^) = 6.5). Hence an additional 
contribution of the gas content from NGC 5954 (Hubble type: Sc) might explain the gas 
surplus since the typical HI mass content for a Sc galaxies is Log^Mffj/D'^) = 6.9. Thus 
we conclude that the ongoing merger already has changed the HI appearance of NGC 5953 
and/or NGC 5954. 



3.1.2. HI Environment and Tidal Interactions 

In order to search for the presence of HI companions we used widefield maps (described 
in Q with a field of view (FOV) of ~ 60' which are sensitive to companions with HI masses 
of > IO^Mq ■ (lower detection limit near the center of the image) and > IO^Mq ■ (at 
radial distances of 30'). Note that the sensitivity scales with the galaxy distance squared, 
indicated as d in units of [lOMpc]. We found that 5 galaxies are surrounded by at least 
one Hl-detected companion with a range in HI mass of (0.01 - 2.8) 10^ Mq. In addition, 
potential optical companions were identified by using NED within our FOV (equivalent to 
a radial distance of 30') and a velocity difference of ±500 km s~^ from the host galaxy (~ 2 
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times the HI velocity bandwith). Optical companions (within these constrains) were found 
for ~56% of our galaxies. All cataloged optical companions (NED, within the constraints 
described above), their systemic velocity difference to the host galaxy and their distance from 
the host galaxy center are listed in Tab. |4| We assume that the relative velocity between 
galaxy and companion is inversely proportional to the disruption time of the galaxy disk. 
For companions with associated HI gas detection (10 companions in total, = 38% of the 
optical companions) also the derived Hl-flux, Hl-mass and their velocity width are provided. 

To estimate the total tidal strength of companions on each galaxy, we used the formalism 



developed by Dahari (1984). The strength of a tidal force affecting a galaxy with the mass 
Mg and the diameter Dq is primarily determined by the mass of the companion and the 
distance Rqc to the power of three between the companion and the galaxy. We approximated 
the absolute distance by the projected distance and the mass of the companion by using the 
relation to its size, Mc a Dl with 7 = 1.5 (|Dahari|[l984| [Verley et aLl|2007D. The tidal 



strength Qgc is defined as the ratio between the tidal force and binding force: 



Ftidai R%c WDgD 



Qgc = — oc oc -3 {i) 

^ bind fyr ^GC 

The sum of the tidal interaction strengths (listed in Tab. |4]) created by all the companions 
in the field is then calculated as Q = logi^^^QGc^- A value of Q < indicates that the 
tidal forces affecting the primary galaxy are smaller than the internal forces, and vice-versa. 

Disturbances were characterized by using the divergence of the rotation curves between 
the approaching and receding side (see for example NGC 4321 in Fig. [?]) as well as obvious 
asymmetries in the HI intensity maps (see for example NGC 1961 in Fig. |3]). A divergence of 
at least So" (statistical error, roughly corresponding to 60 km s^^) of the velocity difference 



\3 



(see ^3.3) in the outer disk between approaching and receding side was required (apart from 
deviations due to streaming motions along spiral arms). From these requirements we found 
that 7 galaxies have a disturbed outer HI disk (the individual galaxies are indicated in 
Tab. [3]). The origin of the disturbances (e.g. tidal forces, ram pressure) and the influence of 



the companions onto the HI disks are discussed in ^4.1 



3.2. Comparison of Stellar and Gaseous Distribution 

In order to compare the stellar distribution with the HI gas distribution, radial profiles 
were obtained by ellipse fitting using the task ELLINT within GIPSY. The radial profiles 
of the mean density and the cumulative sum of the optical light (from DSS2-blue images) 
and HI intensity maps were calculated. The widths of the radii were set to the HI angular 
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resolution which is significantly lower than the typical angular resolution of ~ 2" of the 
DSS2-blue images. In Fig. [i] the mean radial HI and optical profiles are presented (top 
panel) as well as the cumulative distribution of the HI and optical emission (bottom panel). 

To allow for comparison of the spatial extent of the HI gas distribution with the stellar 
component, the ratio between the radius of the HI disk (Rhi) and the radius of the stellar 
disk (Ropt) is calculated. The optical radius Ropt was set to the Holmberg radius defined 
as the parameter log Db25 Q (Hyperleda). In order to obtain an equivalent radius Rhi for 
the HI disk, we measured the radius at a fixed column density of 5.0 x 10^^ cm~^ for all 
our galaxies. The radius of the optical disk (Ropt) and the HI disk (Rhi) are indicated as a 
dashed and a solid line in Fig. |4| respectively. The results are listed in Tab. |3]for all galaxies 
and the mean ratio of is ~ 1.7. In addition the physical radii of the HI disks to this 

limiting column density are given in Tab. [3] with a mean absolute HI radius of 20.6 kpc. 



3.3. HI kinematics 

3.3.1. Kinematic Parameters and Rotation Curve 

The HI morphology provides detailed information about the spatial distribution of the 
gas, but only in combination with their kinematics can the complete dynamical phase space 
of the atomic gas be analyzed. Therefore we have performed a detailed kinematic analysis 
for all galaxies, which includes velocity fields, rotation curves, and residual velocity fields. 
Fig|6] shows (from the left to the right) the intensity-weighted HI mean velocity fields, the 
HI intensity (color) overlaid with the velocity field (contours), the HI velocity dispersion and 
the residual velocity field derived from a smooth fit to the rotation curve. 

The following kinematic parameters (listed in Tab. |5]) have been derived from the ob- 
served velocity field (derived from the robust weighted cubes): 

• The dynamical center and its offset from the optical center (taken from Hyperleda). 

• The systemic velocity (Heliocentric radial velocity) Vgys in km s~^. 

• The position angle (PA) in degrees, defined as the angle between the north-direction on 
the sky and the major axis of the receding half of the galaxy in anti-clockwise direction. 

• The inclination {i) in degrees. 



Length of the projected major axis at the isophotal level 25 mag/arcsec^ in the B-band 
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These parameters were derived by fitting tilted rings to the velocity field using the task 
ROTCUR within GIPSY. Initial starting parameters have been taken from cataloged data 
(NED or Hyperleda). For the fit we excluded data points within an angle of 20° of the minor 
axis. The widths of the radii were set equivalent to the corresponding angular resolution of 
each galaxy (see Tab. [2]). No radial velocity component was fitted. The kinematic parameters 



have been derived in an iterative way following Begeman (1989). The parameters were 
assumed to be the same at all radii (unless specified otherwise), except for the circular 
velocity. Therefore we weighted the obtained value in each ring with its standard deviations, 
in order to obtain the mean parameter. The iterative process of deriving the kinematic 
parameters is as follows: we began by fitting both sides of the velocity field in order to 
determine the systemic velocity (with PA, i held fixed, everything else was free to vary), 
then the kinematic center (with a fixed PA, i, Vgys) and checked again the systemic velocity 
(with a fixed PA, i and center). Next, the derived center, the inclination and the systemic 
velocity were held fixed in order to determine the position angle. The position angle was well 
constrained by the model fits to within ±(l-5)°. The center, the position angle and systemic 
velocity were held fixed in order to obtain the inclination angle. In case of insufficient 
precision (statistical error of the inclination fit was > 8°) for the derived inclination (for 
NGC 2782, NGC 4569, NGC 6574) we applied an ellipse fit to the optical image. The 
inclination derived from the optical axial ratio measured at the Holmberg radius was then 
used to fit the rotation velocity. Finally, the rotation velocities were obtained with fixed 
inclination, position angle, center and systemic velocity. We used for NGC 4826 two different 
PA (inner part: R=0 to 64"; outer part: R=64 to 700") to account for its counter-rotating 
HI disk. For NGC 3718, the inclination was required to increase from z=48° to 60° to model 
its warped HI disk. The derived PA was constant over the disk within ±5°. The derived 
rotation curves for each galaxy are presented in Fig. [Tjfor a) the entire disk (black line), b) 
the approaching side of the disk only (dark grey dots), and c) the receding side of the disk 
only (light grey dots). 

The derived kinematic parameters and rotation curves have been verified by using IN- 
SPECTOR in GIPSY. This task allows us to overlay the circular HI velocities on position- 
velocity (pv) diagrams in order to verify the tilted ring parameters on the 3-dimensional 
datacube. Comparison between the derived rotation curves and corresponding pv slices 
from the data cubes (center and right panel of Fig. [T]) show the effect of beam smearing 
(particularly in the inner 20") and some small misalignments (~ spatial and spectral res- 
olution) between our modeled parameters and observed kinematics (e.g. center, systemic 
velocity). Otherwise, the correspondence between model and data is good in all cases (see 
next subsection). 
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3.3.2. Velocity Fields and Velocity Dispersion 

2-D model velocity fields have been created from the derived interpolated rotation curves 
using the program VELFI within GIPSY. These model velocity fields were subtracted from 
the observed velocity fields to obtain residual velocity fields. The typical range of the residual 
velocities is ± 25 km s~^. Furthermore, the residual velocity fields (Fig. ^ can also be 



utilized to identify systematic misalignments of the derived kinematic parameters (Warner 



et al. 1973). However, no obvious misalignments were identified. 



Most of the galaxies in our sample exhibit fairly regular velocity fields that are dominated 
by circular motion. The three exceptional cases are described in the following: NGC 2782 



exhibits a large gaseous tidal tail (see ^3.1) with a velocity gradient of ~100 km s ^ over 



■^65kpc (360") across that tail. The residual velocity field shows that the HI kinematics in 



the inner HI disk (which coincides with the stellar disk) are fairly different. Smith (1991 



1994) using HI data at ~9" resolution found in addition 1) a shorter HI structure extending 
toward the east (with G.G-IO^M©) and 2) that the HI gas in the inner disk is counterrotation 
with respect to the tidal material, suggesting a merging of two galaxies of unequal masses. 
Due to the resolution of 34" of our HI data, this disk is basically unresolved and the derived 
rotation curve reflects partially the kinematic properties of the tidal tails. 

NGC 5953 exhibits significant streaming motions, inconsistent with global rotation, with 
a velocity gradient of ~60 km s~^ over 100" (~17kpc). Also NGC 1961 shows streaming 
motions in the northern part of its HI distribution (velocity gradient ~400 km s~^ over 
180", ~47kpc), in addition to a rotating disk which is coinciding with the optical one. The 
residual velocity fields revealed for some galaxies of our sample (e.g. NGC 5248, inner disk 
of NGC 6951) indications for non-circular motions. Such motions can be caused by a bar 



potential or spiral waves (van der Kruit & Allen 1978). 



The velocity dispersion maps show that the majority of our galaxies have an average 
velocity dispersion between 10 and 30 km s~^. The reason for the larger dispersion along the 
minor axis seen in almost all galaxies is likely caused by an observational effect: the angular 
distance between adjacent velocities becomes smaller along the minor axis, and hence the 
beam covers a larger velocity gradient. For NGC 3147, NGC 6574, NGC 3627 and NGC 4569 
we found slightly higher dispersion with an average of 15-25 km s~^ and maxima of 40-70 
km s-\ Further, the highly-disturbed galaxies (NGC 1961, NGC 2782, NGC 5953) exhibit 
very large velocity dispersion with maxima of 60-100 km s^^. 
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3.4. Gaseous and Dynamical Masses 

The HI fluxes of the galactic disks were measured in the intensity maps of the naturally 
weighted cubes. The intensity maps were corrected for the response of the VLA primary 
beam. In Tab. [6] the HI flux and the corresponding HI mass are listed for each galaxy. 
The average HI mass of our sample is 5.4 x lO^M© with a maximum mass of 46.6xlO^M0 
(NGC 1961). 



To obtain the dynamical mass we used two approaches. A Brandt curve (Brandt 1960) 
was fltted (black curve) to the measured rotation velocities of the whole disk. The Brandt 
curve is an approximation of a solid-body rotation at small radii and a Keplerian-type velocity 
decrease at large radii and can be described by (see equation [26] and [28] of Brandt) 



viR) R 



V 



Br 



J^Br 
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R 
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where v^^ is the maximum rotation velocity at the radius R^^ and n is an index parame- 
ter 

dyn 



terizing the steepness of the curve. From the derived parameters the dynamical mass M^^ 



M|; = 1127(-) R^^v^^D (5) 



can be calculated (equation [29] of Brandt) via 

'3\3/". 
-2. 

with M^'^ in solar mass units, the distance D in [Mpc], the velocity parameter v^^ in [km 
s~^] and the radius parameter R^^ in ["]. For all galaxies the derived dynamical masses 
are listed in Table |6j In addition, we have calculated the dynamical mass via the spherical 
symmetry approximation: 

Mdyn = MdynlMe] = 232407 Rmikpc] v^[{km s~')% (6) 



with the circular velocity v from the rotation curve at the HI radius Rhi (see ^3.2). Note that 
the mass could be 40% lower, if a disk geometry would be assumed instead of a spherical one. 
The dynamical mass obtained by the Brandt flt and the dynamical mass derived at the 
HI radius Mdyn are presented in Tab. [6j The comparison between both methods shows that 
the mass derived by the Brandt flt is on average a factor of two larger as the mass derived 
at the HI radius. This difference might be expected since the Brandt method is based on 
the assumption that the rotation curve reaches a Keplerian fall off quite rapidly, which is 
not exactly the case for most of our galaxies (probably due to dark matter contribution in 
the outer disk regions). Hence we will use in the following calculations of this paper the 
dynamical mass derived from Eq. |6| The ratio of the HI mass to the dynamical mass of a 
galaxy is given by Mni/Mdyn which lies for our sample in a range of 0.003 to 0.046 with an 
average ratio of 0.015. 
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3.5. HI properties as a function of AGN type 



Here, we investigate possible correlations between the observed HI properties of our 



galaxies and their type of nuclear activity. Since our sample, as described in §2.1[ consists 
mainly of Seyfert (Sy) and LINER (L) host galaxies, we will focus on a comparison of these 
two classes. A summary of the average properties for each AGN type is listed in Tab. [7j The 
error calculation is based on the bootstrapping method ( Efron [1979 ) by taking the standard 
deviation of the mean values from 1000 replications. 



3.5.1. Environment and HI morphology 



As already described in §3.1 we found in our sample several galaxies which are sur- 



rounded by companions and/or have disturbed outer disks. No significant statistical relation 
is found between HI environment and AGN type: The percentage of optical companions is 
roughly the same for each AGN-type (Sy: 3/7, L: 4/7), and the number of galaxies with 
Hl-detected companions (Sy: 1/7 ; L: 2/7; HII: 1/2) is also similar given our small number 
statistics. Galaxies with disturbed outer disks seem to be more abundant in LINER hosts 



(5/7) than in Seyfert hosts (1/7). The origin of these disturbances will be discussed in ^4.1 



In order to search for a possible correlation between the HI morphology and the AGN 
type, we excluded the following galaxies from the analysis, as their HI disk morphology is 
highly asymmetric and hence their morphology can not reliably be assigned: NGC 1961 
(large gas extention to NW), NGC 5953 (interacting system NGC 5953/54), NGC 3718 
(warped HI disk with inner disk close to edge-on), NGC 4826 (very concentrated in center, 
counter-rotating HI disks). Thus the number of Seyferts and LINERs is reduced to 5 galaxies 
for each type. Only LINER host galaxies show ringed HI morphology (LINERs: 80%, Sy: 
%)), while gas spiral arms appear as frequently in both types (Sy: 60%, L: 60%). We 



used the bootstrapping method (Efron 1979) to calculate the statistical significance of this 



result: The confidence interval (at percentiles of [5, 95]%) for LINERs is [40, 100]% for ringed 
structure. The standard deviation of the mean values for HI rings is found to be (80±17)%. 
Furthermore, the radii of the HI rings were determined but neither the absolute radius was 
found to be the same for all ringed galaxies (range between 0.9-5.7 kpc), nor the relative 
radius of the ring to the optical disk {Rrmg/Ropt- 0.15-0.6) or to the HI disk {Rrmg/Ropt- 
0.11-0.6) shows any obvious trend. 
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3.5.2. Masses 

In order to search for a possible relation between the mass properties and the nuclear 
activity type, we compared the dynamical mass (within the HI radius Rhi) as well as the 
HI mass as a function of the AGN type. The mean dynamical mass is slightly larger for 
LINERs (S.S-IO^^Mq) than for Seyfert galaxies {Mdyn is 2.4-lO"M0) but similar within a 
range of [0.45-5.7] -IO^Mq (Sy) and [O.4-12.3]-1O"M0 (L). For the HI gas content, we found 
for Seyfert galaxies a smaller mean HI mass (2.5-lO^M0) than for LINERs (S.G-lO^Mo) with 
a range of [0.2-9.0] -lO^M© and [3.8-46.6] -lO^M©, respectively. In order to study the relative 
gas content we calculated the ratio of MuilMdyn for each galaxy and averaged it for a given 
AGN type. The mean ratio of Mni/Mdyn is 0.009 for Seyfert and 0.013 for LINER host 
galaxies. Interestingly, the two non-AGN galaxies in our sample have a larger HI gas content 
of 0.046 and 0.028. The mean Hubble type is roughly the same for our Seyfert (t=2.8), 
LINER (t=2.7), and non-AGN (t=2.6) sample. 

A comparison between the relative gas content {MH1/D25) versus Hubble type for our 



samples (Sy, L, non-AGN) and typical values for a larger sample (Bettoni et al. 2003) is 



presented in Fig|5j On average the relative gas content of the sample of Bettoni et al. seems 



to be larger compared to our sample. Note that Bettoni et al. (2003) excluded galaxies 
cataloged as having distorted morphology and/or any signature of peculiar kinematics (such 
as polar rings, counterrotating disks or other kinematically decoupled components) and thus 
their sample does not contain galaxies similar to ours. 



3. 5. 3. HI versus optical extent 
In order to compare the relative gas disk size between different AGN types, we averaged 



the ratio Rni/Ropt of the HI radius Rhi to the optical radius Ropt (as described in ^3.2) for 
each AGN type. This analysis shows that the gas disk in Seyfert galaxies is slightly larger 
with an average ratio of 1.9 for Seyferts than for LINERs (1.4). However, the difference 
is within the statistical deviation (± ~0.6). If we neglect galaxies which are part of a 
galaxy group or cluster (NGC3627, NGC4569, NGC4579), the ratios increase to 2.4 and 
1.5, respectively (which corresponds to a difference of 0.9 ± 0.6). Thus, there is a hint that 
Seyfert galaxies might eventually have more extended HI gas disks than LINER galaxies. 
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3.6. Synopsis of results 

The atomic gas distribution and kinematics for a sample of 16 galaxies have been pre- 
sented. There is a large variety of HI distribution in the sample (spiral-, ringed- and centrally 
concentrated geometries) and some galaxies show strongly disturbed HI disks. Most of our 
galaxies exhibit fairly regular HI gas distribution, exceptions are discussed below. The analy- 
sis of the large environment revealed that 5 galaxies have Hl-rich companions and for several 
galaxies potential optical companions are present. The distributions of the gaseous and stel- 
lar component are compared, kinematic properties of the HI gas are analyzed (i.e. rotation 
curves, residual velocity maps), and the gaseous and dynamical masses are derived. Further- 
more, a comparison of HI properties and AGN type suggest correlations with HI morphology 
and host disturbance: LINER host galaxies seem to have a larger fraction of disturbed HI 
disks, a significant higher percentage of HI ring geometry, and tend to possess a smaller HI 
extent (relative to the optical radius) than Seyfert host galaxies. No correlation between the 
presence of companions and the AGN type was found. 



Discussion 



The relationship between large scale environment and nuclear activity has been long 
debated in the literature. Several studies have found indications for correlations between 



the environment and nuclear activity ( Storchi-Bergmann et al. 2001 Chatzichristou 2002 



Marquez et al. 2003). In particular it was suggested that interacting galaxies or galaxies 



with companions exhibit a significant excess of nuclear activity compared to isolated galaxies 



^Dahari 


1984 


Keel et al. 


1985 



seen by other studies (Virani et al. 2000 Schmitt 2001 Fuentes- Williams & Stocke 1988 



MacKenty 1989; Laurikainen & Salo 1995). Hence the issue of possible relations between 
the environment and nuclear activity appears to be still controversial. Furthermore, results 
from Keel et al. (1985) suggested that nuclear phenomena might likely be triggered by a 
tidally induced inflow of gas from the disk to the nuclear regions, rather than gas transfer 
between the interacting galaxies themselves. Keel (1996) also found that Seyfert galaxies 
in pairs actually display smaller kinematic disturbances than non-Seyfert galaxies in pairs, 
which is obviously in disagreement with the hypothesis that tidal interactions are necessary 
for the transport of angular momentum and the fueling of the SMBH. Since most of these 
studies are based on optical/IR imaging, they are in principle less sensitive to distortions 
than our study of the atomic gas that reacts most readily to tidal disturbances. It should be 
pointed out that a detailed study of HI gas properties for active versus non-active galaxies 
(Mundell et al., in prep.) is under- way. 
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In the context of different AGN types, the analysis of a sample of 451 active galaxies 



(Sy, LINER, Transition, HII, and absorption-line galaxies) from the Palomar survey (Schmitt 



2001) showed no correlation between AGN-type and the percentage of galaxies with nearby 



companions after taking morphological differences of the host galaxies into account. We 
also see no evidence for a correlation between the fraction of companions and the AGN 
type present (Seyfert, LINER galaxies), neither from our HI study nor for cataloged optical 
companions listed in NED (see [3.1). Note that our sample has a limited number of LINERs 
(7 galaxies) and Seyfert galaxies (7) and hence the existence of possible relations can not be 
completely excluded. 



4.1. HI environment: tidal forces and their correlation with disturbed HI 

disks 



Companion galaxies can possibly disturb the HI gas in a galactic disk via tidal forces 



and hence affect the fueling of the center with gas (e.g. Keel et al. 1985). The sensitivity 



and endurance of HI to trace the strength and prevalence of tidal interactions among Seyfert 



galaxies is discussed in detail by Greene et al. (2004). In our analysis we found companions 



for about half of our sample. In almost all cases the produced tidal forces are smaller than 



the binding forces of the affected host galaxy (indicated as Q < 0, see ^3.1 and Tab. E 



Only the system NGC5953/54 exhibits signs for very strong gravitational interaction as 



mentioned in ^3.1 Therefore, most of the disturbances in our sample (6/7 galaxies) can not 
be explained simply by tidal forces presently at work. The most probable explanations for 
these disturbances are described in the following, listed in decreasing relevance: 



Interaction with a companion now further or far away. If we assume a relaxation 
timescale of the disturbance of 3-10*yr (typical time for one rotation of a galaxy) and 
a maximum fly-by velocity of ~ 500 km s~^, the involved companion is expected to 
be found within a radial distance of ~ 150 kpc from the disturbed galaxy. Since 71% 
of the disturbed disks show companions in a reasonable distance (projected distance 
is less than 150 kpc), tidal interactions in the past seem to be primarily responsible of 
the disturbances identified in the HI disks. 



Ram pressure stripping (Cayatte et al. 1990 VoUmer et al. 2000 2001). In particular 



for NGC 1961 stripping by intergalactic material was suggested (Shostak et al. 1982). 



But since none of our galaxies with disturbed HI disk lie in a massive group or cluster 
environment, gas stripping due to ram pressure is not very likely to explain the presence 
of disturbances in our sample where no nearby companion is found. Furthermore, the 
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outer HI disk, where disturbances are usually seen, has been removed in clusters due to 
ram pressure stripping, indicated also by a small HI radius (e.g. NGC 4569, NGC 4579 



as part of the Virgo cluster; see ^4.2) 



Minor merging, whereby the companion has now fully merged and has left no optical 
trace. This might be the case for NGC 4736 and NGC 2782 where no companions are 
found in a reasonable distance for tidal interaction. 



Large gas accretion from cosmic filaments: Asymmetries in the gas accretion may 



cause disturbances (for effects of gas accretion on spiral disk dynamics see Bournaud 



& Combes 2002). 



4.2. HI morphology and comparison with the stellar distribution 



The comparison of the radial density profiles between the HI gas distribution and the 
stellar distribution revealed significant deviations as expected: 1.) The extent of the HI 
disk is on average 1.7 times the optical radius, indicated by the Holmberg radius. Only 
NGC 3627, NGC 4569 and NGC 4579 show a smaller HI radius than optical radius. This 
can be explained by the fact that they are all members of interacting groups and/or by their 
rapid motion through an intracluster medium. NGC 3627 is part of the Leo Triplet group 
and the past encounter with NGC 3628 could explain the spatial conincidence of both the 
stars and the gas ( |Zhang et al.|l993| ). The truncated disks of NGC 4569 (and NGC 4579) are 
most probably a signature of strong ram pressure stripping in the past by the intracluster 
medium which pervades the Virgo Cluster (Cayatte et al. 1994). However, most of the 



galaxies in our sample show a larger HI disk than their optical one. 2.) The radial density 
profiles exhibit for most of our galaxies a deficiency of HI gas in the inner part of the galaxy 
disk. This is in general explained by the phase transition from atomic to molecular gas in 



neutral ISM (Young & Scoville 1991). 



4.3. Correlations between HI gas properties and AGN activity type 



As described in ^3.5 our analysis revealed that the number of galaxies with disturbed HI 
disks is higher for LINER galaxies than for Seyfert galaxies. But since several mechanisms 
which can not easily be distinguished can cause these asymmetries (see ^4.1 ), it is not possible 
to draw any strong conclusions explaining the higher fraction for LINERs. 



Our study of the HI morphology revealed a significantly higher percentage of galaxies 
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with a HI gas ring for LINER than for Seyfert galaxies. Interestingly, only the study of 
the Extended 12/im Galaxy Sample ( [Hunt fc Malkan 1999) indicated a prevalence of stellar 



rings in active galaxies (Seyfert and LINERs), where LINERs have elevated rates of inner 
rings, while the Seyfert host galaxies have outer ring fractions several times those in normal 
galaxies. However, we found no HI gas rings in Seyfert host galaxies. Note that stellar rings 
are not preferentially found in Seyfert, LINERs, or non-AGNs for our sample by using the 
optical classification from RC3 listed in NED (see classification in Tab. [l] indicated as R or 
r). To summarize, HI rings are more often found in LINERs while no strong correlation with 
stellar rings is present. 

One possible explanation for an abundance of HI gas rings in LINERs may be a common 
evolution of the gas distribution in the disk together with the nuclear activity where both 
are subject to the influence of a present bar or previous one which has now dissolved. This 
becomes important since rings are linked observationally and phenomenologically to barred 



galaxy dynamics (for a general review see Buta & Combes 1996), and hence, expected to 



be seen after the bar had enough time to redistribute the gas toward the end of its life- 
time. Thus, a time evolution of AGN types seems to be possible, where Seyfert and LINERs 
represent different phases of the galaxy activity cycle: Seyfert galaxies are the ones where 
the fueling process has just been triggered (e.g. through disturbances and/or bar dynamics) 
while LINERs are the ones where the triggering mechanism has already distributed the gas 
in a more stable new configuration (rings) that does no longer support the massive inflow of 
gas. Interestingly, disturbances, which are assumed to be a sign for a recent trigger of gas 
inflow, are preferentially found in LINERs in our sample. That would suggest that LINERS 
are the first stage of activity, just after the triggering: The gas is still distributed in rings, 
and the fueling of the AGN has just started by the redistribution of the gas. However, as 
most of the disturbances in our sample are probably due to tidal interactions in the past, as 



explained in ^4.1, the currently observed nuclear activity must not be identical to the one 



during the ongoing tidal interaction. 

A general scenario for self-regulated activity in low-luminosity AGNs was developed by 



Garcia-Burillo et al. (2005), in which the onset of nuclear activity is explained as a recurrent 
phase during the typical lifetime of any galaxy. In this scenario the activity in galaxies 
is related to that of bar instabilities, expecting that the active phases are not necessarily 
coincident with the phase where the bar has its maximum strength. Since the infall of gas 



driven by a bar is self-destructive (Bournaud & Combes 2002), i.e. it weakens and destroys 



the bar, the potential returns to axisymmetry and the gas piles up in a stable configuration 
(i.e. in rings at the resonances). At this stage, torques exerted by the gravitational potential 
are negligible and other competing mechanisms, e.g. viscous torques, must transport the 
gas in the center of the galaxy. The periods of Seyfert/LINER activity (each lasting ~ 10^ 
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years) are suggested to appear during different evolutionary stages of a bar episode (typically 
characterized by a lifetime ~ 10^ years), depending on the competition between viscosity 
and gravity torques. The prevalence of HI rings in LINERS, derived in this work, could be 



explained by the scenario proposed by Garci'a-Burillo et al. (2005), where AGN activity is 



linked to the evolutionary state of the bar-induced gas flow. However, to further substantiate 
this link requires deriving the complete gravity torque budget based on the HI distribution 
(see Haan et al. in prep.). 

Regarding the relative HI radius {Rhi / RopUcai) we found a slightly larger HI extent for 
Seyfert galaxies than for LINERs, which increases even more when neglecting galaxies which 
lie in galaxy clusters or groups. In contrast, our study of the large environment, the HI 
gas content (by using the ratio Mhi / Mdyn), and the relative HI mass {Mhj/ Mdyn) revealed 
no significant correlation with the AGN-type. This might be due to the limited number of 
galaxies in our sample. Hence we conclude that a detailed HI study with a larger sample 
may reveal more additional information on the interplay between the gaseous component 
and the AGN type present. 



5. Summary 

We present results on the atomic gas distribution and kinematics of 16 nearby spiral 
galaxies. This sample forms the part of the NUGA survey which studies the neutral (atomic 
and molecular) gas dynamics of nearby low luminosity AGN galaxies. We investigated re- 
lations between the atomic gas properties, the environment and the stellar distribution as 
well as possible relations to the AGN type present (Seyfert, LINER). Several galaxies are 
surrounded by companions with associated HI emission and half of our sample has dynami- 
cally disturbed disks. No evidence for a correlation between the presence of companions and 
Seyfert and LINER galaxies is found, neither from our HI study nor for cataloged optical 
companions. 

The HI morphology of our galaxies revealed rings, spiral arms, and centrally concen- 
trated peaks. The comparison with the AGN type present showed that ring structure is 
significantly more often present in LINER than in Seyfert host galaxies. Since bars are dy- 
namically linked to the presence of gas rings, a time evolution of the AGN activity together 
with the redistribution of the neutral gas seems to be a plausible explanation for this finding. 
Additionally, we found a slightly larger HI extent and less prevalent disturbed HI disks for 
Seyfert hosts than for LINER host galaxies. A detailed study with a larger HI sample may 
reveal more relations with respect to the AGN type and would improve the significance of 
our results. 
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Fig. 1. — AGN classification for our sample following the definition of Kewley et al. (2006). The left panel 



shows the separation between AGN (upper right part) and non-AGN (lower left part) galaxies using the 
optical line ratios of [OIIIJ/H/3 versus [NII]/Ha as described in the text. In the right panel the dividing line 
between Seyfert (top part) and LINERs (bottom part) is indicated by using the [OIIIJ/H/3 versus [OI]/Ha 
line ratios (see text). Note that NGC 6574 is excluded in these diagrams, as no accurate ratios of the lines 
used were available in the literature. We used the classification as Seyfert given by NED. 
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Fig. 2. — Channel maps of the natural weighted data cube for the example NGC 6951. 
Each channel has a velocity width of ~20.8 km s~^. The flux cut-off was set at 3a of the 
RMS value (for NGC 6951: 0.0011 Jy/Beam) and the contours are in steps of 3a ■ 2"Jy 
(n=l,2,3...). The channelmaps for all galaxies are available as online data (see Appendix). 
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Fig. 3. — HI intensity maps (upper panel) and the optical emission from DSS2-blue (lower panel with HI 
in contours) of our HI-NUGA sample sorted by AGN type (Seyfert, LINER, HIT). The contours are in n 
steps of 5 • 2" • 10^^ cm^^ of the HI column density (with n = 1,2, 3, ...). 
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Fig. 4. — Comparison between the optical (thin line) and HI (thick line) radial distribution. The mean 
radial HI and optical profile are presented in the top panel and the cumulative profiles of the HI and optical 

emission in the bottom panel for each galaxy. The radius of the HI disk Rhi (thick line) and the radius of 
the stellar disk Rb25 (thin line) are indicated. The values of RB2a are taken from the parameter log(dB2.5) 
found in the literature (Hyperleda) and the radius Rhi for the HI disk corresponds to a measured column 
density of 5.0 x 10^ 



cm-2. 
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Fig. 5. — Comparison of the relative gas content versus Hubble type (Morphology type t) between our 
sample (Sy, LINER, non-AGN) and a larger sample of 1916 galaxies by Bettoni et al (2003 the average 
values are plotted and connected with a solid line). The relative gas content is given as the ratio of the 
atomic gas mass to the square of the diameter in kpc at the isophote of 25 mag arcsec^ Mhi/D^25- 



-37- 




- 39 - 




- 40 - 




- 42 - 




Fig. 6. — Overview of all HI kinematic maps for each galaxy. The intensity- weighted HI mean velocity 
field, the comparison of HI intensity map and velocity field, the velocity dispersion map, and the residual 
velocity field are presented (from the top to the bottom). The contours are set in an interval of 25 km s^^ 
for the isovelocity lines (first and second panel from the top), the central line (thick dark line) represents 
the systemic velocity and white (grey) lines the approaching (receding) velocities. For the dispersion and 
residual field contours are set in 5 km s^^ intervals. The residual velocity field was created by subtracting 
a model velocity field, derived from a smooth fit to the rotation curve, from the observed velocity field (see 
text for details). The white ellipse in the residual field (bottom panel) indicates the maximum radius of the 
model velocity field. 
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Fig. 7. — Presentation of the rotation curve and pv-diagrams for each galaxy. The rotation velocities have 
been plotted in the left panel for a) the whole disk (black dots) , b) the approaching side of the disk only 
(dark grey dots) , and c) the receding side of the disk only (light grey dots) . A Brandt curve ( Brandt|1960 ) is 
fitted (thin black line) as well as a spline interpolation (thick black line) to the derived rotation velocities of 
the whole disk. The center and right panel show the pv diagram of the major and minor axis, respectively. 
The contours are in n steps of the column density at 1.8 • 2" • 10^^ cm~^ (n = 1, 2, 3, ...). Note that the values 
for NGC 2782 are hkely incorrect, as there is a short tidal tail to the east which is barely resolved at our 
resolution (see Smith|[T994 ) . 
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Table 1: Sample Overview 



Name 


RA 


DEC 


Hubble 


Bar 


AGN 


AGN 


Vhel 


Dist 


Res. 




(J2000) 


(J2000) 


Type 




HFS 


Kewley 


[km/s] 


[Mpc] 


[kpc] 


NGC 3147 


10 16 53.65 


+73 24 02.7 


S(rs)bc 




S2 


S 


2820 


40.9 


5.0 


NGC 3627 


11 20 15.03 


+12 59 29.6 


SB(s)b 


X 


L/S2 


S 


727 


6.6 


0.5 


NGC 4569 


12 36 49.80 


+13 09 46.3 


SB(rs)ab 


X 


T2 


S 


-235 


16.8 


1.6 


NGC 4826 


12 56 43.69 


+21 40 57.5 


(R)S(rs)ab 




T2 


S 


408 


4.1 


0.3 


NGC 5953 


15 34 32.39 


+15 11 37.7 


SO-a 




S2 


S 


1965 


36 


3.3 


NGC 6574 


18 11 51.23 


+14 58 54.4 


SB(rs)bc 


X 


S 


S 


2282 


38 


3.3 


NGC 6951 


20 37 14.09 


+66 06 20.3 


SB(rs)bc 


X 


S2 


S 


1424 


24.1 


2.4 


NGC 1961 


05 42 04.80 


+69 22 43.3 


SB(rs)bc 


X 


L2 


L 


3934 


53.9 


4.4 


NGC 3368 


10 46 45.74 


+11 49 11.8 


SB(rs)ab 


X 


L2 


L 


897 


8.1 


0.8 


NGC 3718 


11 32 34.85 


+53 04 04.5 


SB(s)a 


X 


LI. 9 


L 


994 


17.0 


1.1 


NGC 4321 


12 22 54.90 


+15 49 20.6 


SB(s)bc 


X 


T2 


L 


1571 


16.8 


1.4 


NGC 4579 


12 37 43.52 


+11 49 05.5 


SB(rs)b 


X 


S1.9/L1.9 


L 


1519 


37 


1.8 


NGC 4736 


12 50 53.06 


+41 07 13.7 


(R)S(r)ab 


X 


L2 


L 


308 


4.3 


0.4 


NGC 7217 


22 07 52.38 


+31 21 33.4 


(R)S(r)ab 




L2 


L 


952 


16.0 


1.4 


NGC 2782 


09 14 05.11 


+40 06 49.2 


SB(rs)a 


X 


H 


H 


2562 


37.3 


6.2 


NGC 5248 


13 37 32.07 


+08 53 06.2 


(R)SB(rs)bc 


X 


H 


H 


1153 


15 


1.4 



Note. — Summary of the properties of the complete HI-NUGA sample. Listed are only parameters from 
LEDA, NED and the resolution of the VLA observations for all galaxies (beam). The AGN classification 
listed in column (6) is taken from Ho, Fillipenko & Sargent (1997): S - Seyfert, L - LINER, T - transition 
object, H - HII galaxy and NED. The galaxies are divided into separated classes as described in j |2.1| 
Seyfert galaxies(top part), LINER galaxies (mid part) and HII/Transition objects (bottom part), additionally 
indicated in column (7). The velocities listed in column (8) are the assumed systemic velocities. In column 
(10) the spatial resolution of our HI images is listed. 
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Table 2: Observational Setup 



JName 


nat. 


weight 


rob. weight 


KMb nat. 


KMb rob. 


TTITT 7T7T 

b WZl. 




Beam ["] 


Beam ["] 


[mjy/beam] 


[mJy/beam] 


[km / s] 




45.4 


X 


43.6 


29.2 X 


22.0 


0.44 


0.59 


430 




22.5 


X 


19.2 


16.8 X 


15.2 


U.o4 


U.Do 


4oi 


NGC 4569 


41.8 


X 


36.0 


21.2 X 


19.2 


0.48 


0.63 


396 


NGC 4826 


25.6 


X 


21.6 


16.3 X 


14.8 


0.32 


0.28 


372 


NGC 5953 


42.1 


X 


32.5 


22.9 X 


18.3 


0.31 


0.41 


349 


NGC 6574 


30.5 


X 


28.0 


20.2 X 


18.3 


0.31 


0.41 


340 


NGC 6951 


39.6 


X 


34.9 


22.2 X 


18.6 


0.35 


0.45 


338 


NGC 1961 


24.6 


X 


20.1 


18.7 X 


15.0 


0.34 


0.49 


767 


NGC 3368 


44.3 


X 


35.6 


23.0 X 


19.9 


0.28 


0.39 


399 


NGC 3718 


17.4 


X 


15.9 


13.9 X 


13.5 


0.47 


0.84 


465 


NGC 4321 


31.3 


X 


25.6 


18.2 X 


15.9 


0.24 


0.30 


291 


NGC 4579 


43.5 


X 


38.6 


23.1 X 


21.6 


0.33 


0.51 


374 


NGC 4736 


29.2 


X 


28.3 


18.9 X 


15.5 


0.32 


0.28 


284 


NGC 7217 


33.4 


X 


29.6 


18.4 X 


17.0 


0.34 


0.41 


326 


NGC 2782 


52.6 X 48.6 


35.1 X 


33.0 


0.25 


0.37 


246 


NGC 5248 


35.5 


X 


28.5 


19.9 X 


17.5 


0.39 


0.27 


311 



Note. — Overview of the observational properties (VLA) beam size and RMS for natural (nat) and robust 
(rob) weighting. In addition the Full Widths at Zero Intensity (FWZI) of the velocities are listed for all 
galaxies. The velocity resolution is ~20.8 km s~^ (natural) and ~5.2 km s""'^ (robust). Only NGC 2782 has 
a velocity resolution of -^10.4 km for robust weighting. 
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Table 3: HI and Optical Morphology 



Name 


Morphology 


Rhi 


Ropt 


Rhi/ Ropt 


Disturbed 


Companion 




Type 


n 


n 




Disk 


HI detected 


NGC3147 


sp 


200 


122.1 


1.64 






NGC3627 


sp 


193 


267.3 


0.72 






NGC4569 


c 


134 


286.5 


0.47 






NGC4826 


sp/c 


578 


300 


1.93 






NGC5953 


c 


153 


36.9 


4.15 


X 


X 


NGC6574 


c 


71 


45.3 


1.57 






NGC6951 


sp 


255 


97.2 


2.62 






NGC1961 


sp 


130 


125.1 


1.04 


X 


X 


NGC3368 


r/sp 


393 


227.4 


1.73 


X 


X 


NGC3718 


warped 


460 


161 


2.86 


X 


X 


NGC4321 


sp 


237 


228 


1.04 


X 




NGC4579 


r 


147 


168.6 


0.87 






NGC4736 


r/sp 


427 


300 


1.42 


X 




NGC7217 


r 


121 


119.4 


1.01 






NGC2782 


c 


120 


101.7 


1.18 


X 




NGC5248 


sp 


334 


143.7 


2.32 




X 



Note. — Overview of the HI morphology compared to the size of the stellar disk. The HI morphology 

is classified as ring (r), spiral arms (s) and/or bar (b). The values of Ropt arc taken as dB25/2 (Length of 
the projected major axis at the isophotal level 25 mag/arcscc^ in the B-band ) from LEDA. An equivalent 
radius for the HI disk Rhi is defined at a HI column density of 5.0 x 10^^ cm~^. The fraction of the radius 
of the HI disk {Rhi) to the radius of the stellar disk {Rgt) is given as Rhi /Ropt- In addition disturbed disks 
and companions are marked by an "x" (see text for details). 
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Table 6: Mass Properties 



Source 


Flux 


Mm 


Mni/Mdyn 


v{Rhi) 


Rhi 


Mdyn 


A/fBrandt 




[Jy km/s] 






[km/s] 


[kpc] 






NGC3147 


22.7 


8.96 


0.016 


249 


39.7 


5.7 


10.1 


NGC3627 


32.8 


0.34 


0.007 


185 


6.2 


0.5 


0.9 


NGC4569 


7.5 


0.5 


0.003 


245 


10.9 


1.5 


1.2 


NGC4826 


45.9 


0.18 


0.003 


149 


11.5 


0.6 


2.5 


NGC5953 


6.9 


1.77 


0.004 


276 


24.5 


4.3 


4.9 


NGC6574 


2.93 


0.85 


0.019 


126 


12.0 


0.4 


1.5 


NGC6951 


34.8 


4.77 


0.014 


221 


29.8 


3.4 


1.7 


NGC1961 


68 


46.62 


0.037 


402 


34.0 


12.8 


3.4 


NGC3368 


64.2 


0.99 


0.010 


169 


15.4 


1.0 


2.0 


NGC3718 


115.1 


7.85 


0.016 


234 


37.9 


4.8 


8.6 


NGC4321 


51.1 


3.4 


0.013 


238 


19.3 


2.5 


6.6 


NGC4579 


8.5 


0.57 


0.004 


241 


12.0 


1.6 


3.3 


NGC4736 


87 


0.38 


0.009 


141 


8.9 


0.4 


0.2 


NGC7217 


9.7 


0.59 


0.004 


263 


9.4 


1.5 


2.2 


NGC2782 


14.8 


4.86 


0.046 


145 


21.7 


1.1 


1.2 


NGC5248 


81.3 


4.32 


0.028 


166 


24.3 


1.6 


2.0 



Note. — Overview of the mass properties, including the measured HI flux, the dynamical mass -/Vf^™'"** 
obtained from the Brandt fit and the dynamical mass M^y^ derived from the deprojected circular velocity 
v{Rhi) at the HI radius Rhi- The ratio of the HI mass and the dynamical mass Mdyn of a galaxy is given 
by Mhi /Mdyn. Note that Mdyn for NGC 2782 is likely incorrect, as there is a short tidal tail to the east 
which is barely resolved at our resolution (see Smith|[l994 |. 
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Table 7: AGN Type Relations 



Properties 


Seyfert 


LINER 


non-AGN 


Distance [Mpc] 


22.9 


20.0 


26.2 


Morph. Type t 


2.8 ± 0.5 


2.7 ± 0.4 


2.6 


Disturbed (%) 


14 ± 13 


71 ± 17 


1 


HI detected companions (%) 


14 ± 13 


43 ± 18 


1 


Rhi [kpc] 


19.2 


19.6 


23.0 


Rb25 [kpc] 


12.4 


14.7 


14.4 


Rhi 

Rb25 


1.9 ± 0.43 


1.4 ± 0.25 


1.8 


HI Mass [IO^Mq] 


2.5 ± 1.1 


8.6 ± 5.9 


4.5 


Mhi 

^dyn 


0.009 ± 0.002 


0.013 ± 0.004 


1.8 


Ring (%) 





80 ± 18 





Spiral (%) 


60 


60 


100 



Note. — Overview of the mean HI properties as a function of AGN type, for Seyferts, LINERs and non- 
AGN galaxies. For each AGN type the number of galaxies, the average distance, the percentage of galaxies 
with disturbed disk and HI companions, the mean HI radius Rhi at a HI column density of 5.0 x 10^^ cm^^, 
the mean optical radius Rb25, the ratio Rhi / Rb25i the average HI mass of the disk, the average morphology 
type t, and the percentage of ringed, spiral and concentrated HI structure are listed. In addition, also the 
corresponding standard deviations are specified (using Bootstrapping method, see ^3.51. 



